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Abstract. The dilepton production in elementary pp — > e~^e~ X reactions at 
BEVALAC energies Tiab = 1 -r 5 GeV is investigated. The calculations include 
direct e+e~ decays of the vector mesons p", oj, and </>, Dalitz decays of the tt"- 
, rj-, p-, UJ-, and 0-mesons, and of the baryon resonances A(1232), A'^(1520), ... 
. The subthreshold vector meson production cross sections in pp collisions are 
treated in a way sufficient to avoid double counting with the inclusive vector 
meson production. The vector meson dominance model for the transition form 
factors of the resonance Dalitz decays R — > e'^e~ N is used in an extended form to 
ensure correct asymptotics which are in agreement with the quark counting rules. 
Such a modification gives an unified and consistent description of both R — > N-y 
radiative decays and R — > Np{Lu) meson decays. The effect of multiple pion 
production on the experimental efficiency for the detection of the dilepton pairs is 
studied. We find the dilepton yield in reasonable agreement with the experimental 
data for the set of intermediate energies whereas at the highest energy Tj^i, = 4.88 
GeV the number of dilepton pairs is likely to be overestimated experimentally in 
the mass range M = 300 700 MeV. 
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1. Introduction 

Relativistic heavy-ion collisions present an unique possibility to create nuclear matter 
at high densities and temperatures where the hadron properties become different and a 
phase transition to quark matter with signatures of the deconfinement and restoration 
of chiral symmetry is expected. The change of the nucleon mass in the nuclear matter 
was implemented into the Walecka model 0E| already in 1970's in the framework of 
the effective hadron field theory. Later on this effect was put on the firmer grounds 
on the basis of a partial restoration of the chiral symmetry and finite-density QCD 
sum rules ^ , the Nambu- Jona-Lasinio model [S] , and effective meson field theory 

Dileptons (e+e^, /i^/i^) are the most clear probes of the high density nuclear 
matter. The reason is that the dileptons interact with the matter only by 
electromagnetic forces and can therefore leave the heavy-ion reaction zone essentially 
undistorted by final state interactions. They provide valuable information on the 
in-medium properties of hadrons and, hence, on the state of matter. 

The dilepton spectra from heavy-ion collisions have been measured at two different 
energy scales: by the CERES and HELIOS-3 Collaborations at SPS E] (a few 
hundreds GeV per nucleon) and by the DLS Collaboration at BEVALAC ^01 (a few 
GeV per nucleon). In the CERES and HELIOS-3 experiments and in the BEVALAC 
experiment the production of dileptons with invariant masses between 300 700 MeV 
is found to be enhanced as compared to estimates based on the theoretically known 
dilepton sources when in-medium modifications of hadron properties are neglected. 

The data on the total photoabsorption cross section on heavy nuclei [11, give an 
evidence for the broadening of nucleon resonances in the nuclear medium |12:. The 
physics behind is the same as in the collision broadening of the atomic spectral lines 
in hot and dense gases, discussed by Weisskopf in the early 1930's (13j (for the present 
status of this field, see, e.g., Jl]). This rather general effect which takes its origin 
from the atomic spectroscopy should also lead to a broadening of the p-mesons in 
heavy-ion collisions. 

The low-energy dilepton excess can be explained by a reduction of the p-meson 
mass in a dense medium [151 1161 117) . The in- medium modification of the p- meson 
spectral function jl8[ [T^ leading to a broadening of the p-meson seems also to be 
sufficient to account for the CERES and HELIOS-3 data [Tt] . 

In the DLS experiment a different temperature and density regime is probed. The 
enhancement of the dilepton spectra due to the reduction of the p-meson mass and 
the p-meson broadening is not sufficient to bring theoretical estimates in coincidence 
with the available experimental data '20". The in- medium p- meson scenarios that 
successfully explain the dilepton yield at SPS energies fail for the DLS data. This 
phenomenon was called the " DLS puzzle" . It worthwhile to notice that the final data 
from the DLS Collaboration have changed by about a factor of 5 — 7 as compared to 
the initially reported results. The future HADES experiment at GSI will study the 
dilepton spectra at the same energy range in greater details |21| . 

A possibility to clarify the origin of the DLS puzzle has appeared since data 
from elementary pp (pd) collisions at T = 1 -f- 5 GeV (T is the kinetic energy of the 
incident proton in the laboratory frame) became available from the DLS Collaboration 
|22| . The elementary cross sections enter as an input into the transport simulations 
of heavy-ion collisions, so their better understanding is of great value. 

The dilepton spectra in the pp collisions at T = 1-^5 GeV have been calculated in 
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refs. |23[ 1^125) . In ref.|23| the agreement achieved with the DLS data is generaUy good 
at low energies where the subthreshold production of nucleon resonances is important. 
When the energy increases and the inclusive production becomes dominant, the 
dilepton yield is underestimated at the same mass range 300-^700 MeV as in the heavy- 
ion collisions. A signature for this effect exists also in the calculations of refs. |24ll^ . 
This can be interpreted to mean that the studies |23l 1241 12 5| revealed, apparently, the 
reoccurrence of the DLS puzzle on the elementary level of the nucleon collisions. They 
leave, therefore, a doubt on the quality of the experimental data and/or the reliability 
of the accepted theoretical schemes. 

This paper is devoted to a further going theoretical analysis of the elementary 
dilepton production cross sections. 

In the next Sect., the production mechanisms are critically revisited. The 
subthreshold production cross sections for the vector mesons are treated such that 
no double counting appears with the inclusive processes. The effect of multiple pion 
production on the experimental detection efficiency for the dilepton pairs is also 
studied. We demonstrate that the detector efhciency is sensitive to the number of 
pions produced and propose a simple model to account for the multiple pion production 
effects. 

In the Dalitz decays of the nucleon resonances, R Ne^e~ , the Vector Meson 
Dominance (VMD) model is usually applied for the description of the resonance 
transition form factors. However, the naive VMD which takes the R Np data 
as an input systematically overestimates the radiative decay branchings. In Sect. 3, 
the VMD model is extended to ensure the correct asymptotic behavior of the transition 
form factors in agreement with the quark counting rules. Such a modification is found 
to be sufficient to achieve an unified description of both, R radiative decays 

and R N p{u>) meson decays. Our estimates of the subthreshold cross sections rely 
therefore on the two essentially different sets of the experimental data, R N'j and 
R Np. The numerical results are discussed in Sect. 4. We found that the above 
improvements do not eliminate the discrepancy with the DLS data at T = 4.88 GeV. 
Moreover, the results for the lowest energy, T — 1.04 GeV, also require an additional 
study from the experimental and/or theoretical side. 



2. pp ^e+e X reaction 



The dilepton production in nucleon collisions goes through the production of virtual 
photons which decay subsequently into e+e~ pairs. According to the VMD model, 
the virtual photons are coupled to vector mesons V = p^, to, and (j)- The dilepton 
production can therefore be calculated using the inclusive vector meson production 
cross sections: 

Here, s is the square of the invariant mass of two colliding protons, M the invariant 
mass of the dilepton pair, da{s, AI)pp^^^ /dM'^ is the differential vector meson 
production cross section, and ny is the average number of additional vector mesons V 
in the state X. In the energy range of interest, T = 1 5 GeV, where T is the kinetic 
energy of the proton in the laboratory system, ny = 0. The branching ratio 
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corresponds to the direct V — > e+e^ decays, with T^^^{M) being the total meson decay 
width. 

In order to disentangle the various contributions, we decompose the cross section 
entering into eq.(E3) into pole and background parts: 

da(s, = da{s, MYp^ + da{s, . (3) 

Such a decomposition implies that interference effects between the different sources 
are neglected. The background sources like tt — > 7e+e~ and 77 — + 7e+e^ do 
not interfere due to the kinematical reasons. These two reactions in turn do not 
interfere with the reaction u) 7r°e+e~, since the final states are different. The 
reactions NN — > RN , R — > e+e^X going through the nucleon resonances R with 
different quantum numbers do not interfere with each other either. The sources like 
NN A^A(I232), A(1232) Ne+e" and NN NNp, p e+er do interfere, 
however. The relative phases of the amplitudes describing the different reactions are 
unknown, so the neglection of the interference between all the reactions constitutes a 
reasonable first approximation. 

The distribution over the meson mass M in the pole part of the cross section has 
a Breit-Wigner form corrected to the available phase space for the final state VX. At 
moderate energies, the state X is dominated by two nucleons and pions, so one can 
write 

1 MrL(M)dM2 



do{s,M)V^o{s)p --— — ^ 



''VI 
N 

^ w„C„$3+„(Vs--) (4) 

n=0 



where 



$3+ri(\/s...) = ^{^/s,mN,mN,M,flj,, ...,M7r) (5) 

is the (3 + n)-body phase space of the final state (two nucleons with masses m^r, 
one vector meson V with mass my, and n pions with masses fiTt). The value 
iV^r = [{y/s — 2mN — my) / fJ-Tr] is the maximal number of pions allowed by energy 
conservation, [x] denotes the integer value of x, and the values Wn are the probabilities 
for the production of n pions, with 

J2w^ = l. (6) 

The normalization factor C„ is given by 

_ riV-s-2,n.-n^.f 1 MTYAM)dM^ 

where /^o is the physical threshold for vector meson decays (/io = 2/i7r for the p- 
meson). Notice that the cross section vanishes at values M < pQ. However, the 
total width VYot (M) entering into the denominator of the branching ratio at Af < 
Pq vanishes as well, so that the cross section ( is actually finite everywhere above 
the two-electron threshold. 

In the zero- width limit, TYot{M) = 0, eq.Q simplifies to give 

da{s, MYp^ = a{s)Y,^5{M^ - mlr)dM'^. (8) 
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The finite-width effects are important for the p-meson and less important for oj- and 
(/)- mesons. 

As the background we consider states with 7r-mesons originating from the V- 
meson strong decays, that do not contribute to the pole part of the cross section of 
meson production. We assume e.g. that the p- mesons from the reactions 



^0 






{-yn+TT ), 


■q - 




-> 7e^e^ 


(77r+7r"), 






> n^e+e- 


(7r°7r+7r-) 


R- 




-> Ne+e- 


- (Ntt+tt-) 



form a background, which should be added to the pole part of p-mesons production. 
(In the last line only the p-mesons away from the resonance region correspond to 
background.) 

The dilepton invariant mass in all these reactions can be lower than the two-pion 
threshold. Clearly, all processes with M < 2p7r are not accounted for by the pole part 
of the inclusive cross section. Their contributions should be added to the contribution 
of the pole part. With increasing invariant dilepton mass, somewhere in the region 
M w mp, we come to the double counting problem. This part of the spectrum from 
the above reactions must be excluded as outlined below. 

Experimental data on the exclusive cross sections a{s)'^^ with X — rnrNN at 
n > 1 are not available. Here we assume that the probabilities Wn are described by a 
binomial distribution 



n!(A'7r — ny. 

To fix all probabilities it is sufficient to know the ratio between the exclusive vector 
meson production cross section cr(s)p^^ and the inclusive cross section (t(s)p"^. These 
two cross sections are experimentally known |27j . The value is defined as above 
by energy conservation while the value p can be extracted from the relation 



{l-pr-. (10) 



IP 

The pion multiplicity equals 

71^ = ^ nwn = pN^,. (11) 
ji=0 

In the case of the p-meson, the cross section a{s)Y^ determines the pole behavior 
of the total cross section da{s, MY^'^ ^ in the vicinity of the p-meson peak. Like for 
vector mesons, eq.Q, the total cross section da{s^MY^^ ^ for the 27r production 
can be decomposed into p-pole and background parts 

da(s, MY^""'^ = d(T{s, Mfp^ + da{s, MYb""'^- (12) 

The background parts on the right hand sides of ea. H12|l and eq.© do not coincide 
since the 7r"'"7r~ quantum numbers are not necessarily equal to the quantum numbers 
of the p-meson (this is the case at the p-meson peak only and this is why with V = 
the first terms in eq.© and ea. H12|) coincide). Thus the left hand sides of ea. H12|) and 
eq.||21) do not coincide as well. Since the p-meson is always detected via 2-k final states 
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(or via Stt and 2K final states for uj- and 0- mesons, respectively), the inclusive cross 
section for the production of vector mesons which enters into eq.([3) cannot be uniquely 
determined from the experimental data on inclusive production of tt+tt" . Instead, 
one needs a model for the calculation of the background part of the cross section in 
eq.©. A subtle problem of double counting in the total dilepton production cross 
section appears in this way. We propose a quite natural phenomenological solution 
for it. 

The background term da{s, M)q a.t X ^ NN can be saturated, at least 
partially, by considering the production of light mesons: pp —> rjX — > p^jX —> 
Tr~^Tr~jX, pp — > LuX — > p^tt'^X —>■ tt^tt^tt'^X, etc., similarly for the uj- and (/)-mesons. 
The TT+vr^ invariant masses are small here, so these processes contribute to the tt+tt^ 
background. Eq.QJ can now be rewritten as follows 

da{s, Mf^-''' = da(s, Mf^ + da{s, Mfs"''' \x=nn 

+ do{s,MYs''''\x^NN- (13) 
The first term is the same as in Eq.(^), i.e. with the sum running over all vector 
mesons, however, keeping only the pole part of the cross section. In Ea. (|13|) the 
background is divided into contributions from direct decays of intermediate vector 
mesons, which are off-shell and typically below their physical thresholds {X = NN) 
and from decays of intermediate mesons, A4, to multi-particle final states {X ^ NN). 
In the latter case, the cross section for the production of the intermediate meson has 
to be folded over its branching ratio to the final state under consideration: 



dM2 '^^^^ 

M 



da{s,M)<g<''^ . V- /"j 2/. 

^ I ^zZ d^Ji^{l + nM) 

M 

d<j{s,fi)^^' dB{fi,M)^-^^^^" 

dfj."^ dAP ' ^ ' 

The sum runs over the mesons M = tt, rj, p, lu, and 0. Here, is the average 
number of mesons A4 in the state X' . The value p in the last equation describes 
the distribution over the off-shell masses of the mesons. For pseudo-scalar mesons, 
the cross sections due to their small widths are proportional to the delta-function 
d{iJ,^ — rnj^) (cf. Eq.JHl) and the expression reduces to a sum over the on-shell mesons 
decaying to the states e'^eX" with X" ^ Cd, X = X' + X". For vector mesons entering 
the sum of Eq. (|14|l one should use the cross sections JSJl whose pole components are 
well defined. 

The contribution to the background part of the cross section with X = 
NN can be calculated assuming that it results from subthreshold decays of baryon 
resonances R = A(1232), N{1520), ... produced in pp collisions, which decay into 
nucleons and vector mesons, R — s- NV 25 . In terms of the branching ratios for the 
Dalitz decays of the baryon resonances, the cross section can be written as follows 



^M2 \X = NN ^ / "^A^^ J, ,2 



X 



^ dB{p,M)^-^''P-- 

^ dKP 

V 
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Here, /i is the running mass of the baryon resonance R with the cross section 
da{s, fi)PP^P^^ dB{p., M)^^^P"''=^"="P is the differential branching ratio for the Dahtz 
decay R — s- e+e^p through the vector meson V . 

With increasing energy, the vector mesons in Eg. ljlSI) can be produced at their 
physical masses. In such a case, the processes (|15|) contribute to the pole part of the 
cross section da{s^ MY^ ^ . The pole part of inclusive cross sections by definition 
accounts for all possible sources for the appearance of on-shell vector mesons, so that 
a naive extension of the subthreshold cross section to higher energies would result in 
a double counting. To avoid this double counting we should skip the exclusive part of 
vector meson production in pp collision that enters to the total inclusive production 
of vector mesons as already taken into account in Eq. H15() : 

da{s) = da{s)"''\\ - Wo) + dcr(s)""''*'' (16) 

The factor 1 — wo excludes the NNe'^e'^ component of the inclusive cross section 
which is reproduced by the subthreshold mechanizm. This prescription lies on the 
comparison of the experimental data to the exclusive vector meson production cross 
section, that goes through baryon resonances 



^{s) \x=NN ^ 2^ I 



dB{ii,M)^^ ^P 
dJVP 



divr . (17) 



The cross sections reasonably agree (see Figs. ^andlJll. 

The meson multiplicities n» in Eq. 114|l are set equal to zero except for the neutral 
pions. For pions, we assume n^o = injr where n-^ is given by Ea. (|ll|l . The factor i 
appears statistically in the limit n — > cx) from charge conservation which implies that 
in every channel the total number of the neutral-pion pairs tt'^tt^ must be equal to the 
number of the charged-pion 7r+7r~ pairs. 

There are no experimental data on the pion production inclusive cross section 
a{s, ijLTr)PP^^ at r = 1 5 GeV. There exist, however, experimental data on the 
two-pion production cross sections a{s, ^^^p^'^'^^^ |2S|- The parameter p of the pion 
number distribution can therefore be estimated from equation 

^(^)^""^"^ ^ (18) 

(j{sY"pp i-p' ^ ' 

It can further be used to find the total inclusive cross section (t(s, ^i-kYP'^^ ^ using 
the same relation as in Ea. HlU|) . The results for the pion multiplicities produced 
in reactions in accordance with different mesons M. are summarized in Table 1 for 
those values of kinetic energies, T, at which the dilepton cross sections are measured 
by the DLS collaboration. We give there also the cross sections obtained from the 
interpolation and/or extrapolation between the available experimental points. The 
accuracy of these estimates is rarely better than 20%. For the pion, we give an 
estimate for the inclusive cross section. Notice that it is in reasonable agreement with 
the prediction from the UrQMD transport model |23|. 

For mesons M — rj, p, uj, and (p, estimates for the cross sections pp MttNN 
are given as derived from the distribution 

Also the estimates for the exclusive multi-pion production cross are in line with 
prediction from the FRITIOF string fragmentation model '5^ and the calculation of 
ref.jSni- To show this we make the following comparisons: 
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Ec.m. [GeV] 

Figure 1. Exclusive meson production in proton-proton collision 

calculated over contribution of intermediate baryon resonances 1171 compared 
to experimental data. 



First, we compare our simple formulas explicitely to FRITIOF. In refs.[21IS2 an 
attempt was made to extrapolate the FRITIOF model, originally proposed for high 
energies, to lower energies by a modification of the two-body mechanism for inelastic 
hadron-hadron reactions. The authors claimed that they got finally a reasonable 
description of the basic characteristics of the hadron-hadron and nucleus- nucleus 
collisions, within the modified FRITIOF model. 

We compare our model with the modified FRITIOF model of refs.|Hll 02] for 
the reaction np ppir^X at Piab — 3.83 GeV and 5.1 GeV. In Table 2, we show 
4 channels with X = 0, tt", tt+tt", 7r+7r~7r''. The cross sections for the channels 
X — TT^TT^, TT^TT^TT^ arc uot available. We assume 

a{np — s- pptt^tt^tt'^) — a{np — > ppn^ tt^ ) , (19) 

cr(np — !■ ppn^ tt'^ tt'^ n'^) = a{np ppn^ n'^ ) . (20) 

We use the experimental data from the first and the last lines of the Table 2 to fit 
the parameter p of the binomial distribution. The results are shown in the Table 2 
(third lines in the boxes), together with the experimental data (first lines of the boxes, 
numbers with errors) and predictions of the modified FRITIOF model (second lines 
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Table 1. The maximum number of pions, A^Tr, the average pion multipUcities, 
n,r, and the cross sections Q-M-irNN ^ ^MNN ^ g^j^^^ a-^-^ for production of mesons 
A4. = TT^ , r], p, u!, and (p in the proton-proton coUisions for the set of energies 
T = 1.04, 1.27, 1.61, 2.09, and 4.88 GeV, at which the dilepton production cross 
section has been measured by the DLS Collaboration. The last three lines with 
the cross sections contain 3 X 12 = 36 numbers. The 24 numbers are compillation 
of the available experimental data, while the 12 numbers marked with the symbol 
"#" 3^j.g predictions from our binomial formula j^. The experimental cross 
sections are measured at energies different from the BEVALAC energies, so we 
give interpolations from the available experimental points. 



T [GeV] 


1.04 


1.27 


1.61 


1.85 


1.85 


2.09 


M 


ttO 


ttO 


tt" 


tt" 


V 


tt" 




2 


3 


3 


4 


1 


5 




0.09 


0.20 


0.56 


0.74 


0.26 


1.07 




0.4 


0.9 


2.7 


3.4 


4.9 


0.05#) 


^MNN [j^b] 


4.5 


4.2 


3.9 


3.8 


0.13 


3.6 


cr^^ [mb] 


4.9#) 


5.2#) 


7.3#) 


8.5#) 


0.17 


12#) 


T [GeV] 


2.09 


4.88 


4.88 


4.88 


4.88 


4.88 


M 


ri 




V 




LU 






2 


11 


8 


6 


6 


4 




0.55 


1.80 


1.87 


1.67 


1.36 


1.54 




o.ii#) 


5.8 


0.34#) 


0.28#) 


0.35#) 


0.003#) 


^MNN [j^b] 


0.14 


2.7 


0.14 


0.12 


0.20 


0.001 


cr-^^ [mb] 


0.27 


19#) 


1.19 


0.94 


0.94 


0.007 



Table 2. The multipion production cross sections estimated in the FRITIOF 
model (second line) and from the binomial distribution (third line) for the proton 
beam momentum of Piab = 3.83 and 5.1 GeV. Given in the first line are the 
experimental data. 



Plab [GeV] 


3.83 


5.1 


np — > ppn^ 


2.35 ±0.12 
(1.62) 
2.35 


2.13 ±0.11 
(1.75) 
2.13 


np pptt~tt'^ 


1.83 ±0.13 
(3.03) 
2.05 


2.05 ±0.12 
(2.60) 
2.99 


np — > ppTT^ TT^ 


0.31 ±0.04 
(0.41) 
0.38 


0.56 ±0.04 
(0.90) 
0.93 


np ~~> pp-n^ TT^ TT^ TT^ 


0.08 ±0.01 
(0.06) 
0.08 


0.34 ±0.03 
(0.28) 
0.34 
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Figure 2. Exclusive u) meson production in proton-proton collision 
calculated over contribution of intermediate baryon resonances 1171 compared 
to experimental data. 



of the boxes, numbers in the brackets). The experimental data are from ref. . Is 
is seen that the description is reasonable. Notice that a 30% accuracy occur for the 
modified FRITIOF predictions and similarly 30% is also the accuracy of the binomial 
formula predictions. 

Second, the agreement with the work j30| is reasonable. For example as it can 
be seen from Table 2 that our value for cr{np — > ppn^n^) at 5.1 GeV is close to their 
result of 2.85 mb at 5.5 GeV. 

In Table 3, we compare predictions from the binomial distribution for the reaction 
pp — !■ ppn'^X with the experimental data from ref.|!2B . The uncertainty exists in these 
data also: the cross sections for X — 7r°, tt^tt", tt^tt^tt", tt^tt+tt", 7r°7r°7r"7r", and 
TT^TT^TT^Tr"'" are unknown. For a qualitative estimate, we assume like before that the 
cross sections for the unknown channels are equal to the cross sections for the known 
channels, with the same number of the pious. The parameter p of the binomial formula 
is determined at each energy from the channels with one pion and with the maximum 
number of the pious in the final states. 

The agreement is not unreasonable, especially in the view of the naive 
approximations like (|19|l and H20() . The binomial formula does not contradict to the 
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Table 3. The multipion production cross sections in proton- proton collisions 
estimated form the binomial distribution for different momenta of the beam proton 
(first coulomn). The second coulomn shows the final states, the third one shows 
the parameter p of the binomial distribution, the next coulomn gives the estimated 
cross sections in mb, and the last one gives the experimental numbers in mb. 



Plab [GeV] 


final state 


P 




^expt 


5.5 







127 


2.77 


2.77±0.11 


0.0 


ppTT^ TT 





127 


A 0/1 




5.5 







127 


1.327 


1.81 ±0.07 


5.5 







127 


0.516 


0.227 ±0.023 


5.5 







127 


0.088 


0.088 ±0.014 


3.68 


ppir'^ 





134 


2.950 


2.95 ±0.31 


3.68 


ppTT^TT^ 





134 


3.213 


2.72 ±0.13 


3.68 


pp-K^TT^i:^ 





134 


0.750 


0.75 ±0.07 


2.81 


ppir'^ 





095 


3.600 


3.60 ±0.21 


2.81 







095 


1.897 


2.35 ±0.14 


2.81 


ppTT^TT^TT^ 





095 


0.200 


0.20 ±0.03 


2.23 


ppn'^ 





054 


4.060 


4.06 ±0.27 


2.23 


ppTT^TT^ 





054 


0.697 


1.24±0.14 


2.23 


PPTT^TT^TT^ 





054 


0.020 


0.02 ±0.02 



available data and we suppose that it can be used for estimates of the multiple pion 
production effects. Notice that the pion contribution to the dilepton spectrum is quite 
far from the really interesting region about and just below the p-meson peak. 

For the mesonic dilepton decays of mesons M e~^e~X, experimental data 
exist in most cases. The radiative decays Ai — > can further be used to calculate 
the dilepton decays when the experimental data are not available. For details of the 
determination of the various branching ratios see ref . |34| . The inclusive cross sections 
for the meson production and the branching ratios for the mesonic decays to dileptons 
can be combined to estimate the dilepton yield from Ea. H13|l . 

The differential decay branchings dB^^^*^^ ^ {p,, M) / dM"^ are calculated in 
ref. j25| in a non-relativistic approximation for the multipole decays with the emission 
of a massive vector particle. We follow a similar approach here, but consider the 
relativistic case and modify the transition form factors for the nucleon resonances, 
R, which is needed to bring their asymptotic behavior in the correspondence with 
the quark counting rules and to provide an unified description of the photo- and 
electroproduction data and vector meson decays R Np{uj) ,38, . 

3. Transition form factors, quark counting rules, and radiative and vector 
meson decays of nucleon resonances 

The problem is adequately formulated in the non-relativistic approximation for 
radiative and vector meson decays of nucleon resonances. We start with the discussion 
of this case, since it is much simpler. The relativistic treatment, which will finally be 
implemented into the consideration, although is more complicated, is motivated by 
the same physical ideas. 
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Table 4. The coupling constants fuNp derived from the R N p mesonic decays 
are compared to the coupling constants f]ij^p fixed from the radiative R — > N'y 
decays. The numerical values fuNp are taken from ref. with exception of the 
A(1232) resonance for which the theoretical value from 1181 is given and of the 
A''{1440) and A''{1535) resonances where the results of our calculations arc given. 



R 




^1520 


iVl535 


^1650 


^1680 


-^1720 


Ai232 


Ai620 


Ai700 


^1905 


JP 


1 + 
2 


3- 
2 


1 - 
2 


1 - 
2 


5 + 
2 


3 + 
2 


3 + 
2 


1 " 
2 


3- 
2 


5 + 
2 


fRNp 


< 26 


7.0 


< 2.0 


0.9 


6.3 


7.8 


15.3 


2.5 


5.0 


12.2 


fl 

JRNp 


1.3 


3.8 


1.8 


< 0.8 


3.9 


2.2 


10.8 


0.7 


2.7 


2.1 



3.1. Vector meson decays of nucleon resonances in the nonrelativistic approximation 

The description of the resonance decays R A^7*, 7* e'^e~ is usuaUy based on the 
VMD model which provides transition form factors RNj of a monopole form. The 
pole corresponds to the masses of the p- and cj-mesons. This model should give, in 
principle, an unified description of the radiative RNj and the mesonic RNV decays. 
However, a normalization to the radiative branchings (RN^) strongly underestimates 
the mesonic branchings {RNV) as we discuss below. 

The resonance A^(1520) is a case for which both, the A^(1520) Np and 
-/V(1520) widths are known with a relatively high precision: i?(7V(1520) 

Np) = 15 ^ 25%, S(7V(1520) ^ Nj) = 0.46 ~ 0.56 % (pj mode), 0.30 0.53 % {nj 
mode) . The branching ratios of the proton and neutron modes are equal within the 
experimental errors. Let us, for the moment, interpreete this by assumption that the 
radiative mode is dominated by the p-meson. The same conclusion is reasonable for 
other N* resonances: B{N{1U0) iV7) = 0.035-^-0.048 % (p7 mode), 0.009-^-0.032 
% (n7 mode); B(iV(1535) Nj) = 0.15 0.35 % (p7 mode), 0.004 0.29 % {nj 
mode), etc. The A decays proceed exclusively through the p- meson. 

However, now the standard VMD model as it has been used in leads to a 
severe inconsistency: Using the coupling constant fN{i520)Np = 7-0 extracted from the 
mesonic iV(1520) Np decay, the branching ratio for the radiative decay is found to 
be two to three times greater than the experimental value. Analogous overestimations 
are observed almost for all other N and A resonances for which the experimental Np 
and A^7 data are available. Table 4 summarizes the results. 

The standard VMD predicts a 1/t asymptotic behavior for the transition form 
factors. However, quark counting rules require a stronger suppression at high t. It 
is known from the nucleon form factors, the pion form factor, and the ton^ and pirj 
transition form factors that the quark counting rules start to work experimentally 
at moderate t 1 GeV^. One can assume that an appropriate modification of the 
standard VMD which takes the correct asymptotics of the RN^ transition form factors 
into account can provide a more accurate description of the radiative decays of the 
nucleon resonances. 

We propose the following solution of the inconsistency between the RNV 
and RNj decay rates: Let radial excitations of the p-meson, the /9(1450)-meson 
and /5(1700)-meson, for example, interfere with the p- meson in radiative processes. 
However, we know neither the couplings of the p{lA5Q) and /9(1700) to the resonances 
ifRNp', fRNp") nor the coupHngs of the p(1450) and p(1700) to a photon (^p/, (?p"). 
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Thus in the sum 



^ ' mp gp rhf,- rUp, gp, m^, - 

_^ /rjvp" ml„ 1 



P" gp" m^,, 

iMTk with k — p, p' , and p", p' and p" refer to /9(1450)- and 
/c»(1700)-niesons, respectively, the coefficients "gi~ ^"^^ m"» "g"T7 ^'^^ unknown. 

According to the quark counting rules |4()l HT] , for large and negative the form 
factors of the RN^* amplitudes decrease like On the phenomenological level 

we can attribute such a behavior to a cancellation between the p-, p'-, and /^''-mesons. 

The constants and i^EjiL^ElL are then fixed and we obtain 



Jrnp 



m 




M{M^) ^ '-^-^-r,-^, I ^1^^ 1 I 1 . (21) 

The last two factors in Eq.|2U give the desired modification of the p-meson 
contribution to the radiative decays of the baryon resonances, as compared to the 
naive VMD model: 

The mass-dependent correction factor is given by 



Fp{M' 



— 




(23) 



-M2 

The same modification applies to the R decays. The reduction factor in the 

amplitude R — > iV7 equals ^jFpiJXP = 0) = 0.56. It is seen from Table 4 that a 
reduction of about ^ is just what one needs for a consistent description of both, the 
p-meson and the radiative decay of the Af(1520). In all other cases the reduction 
factor also improves the agreement. In the case of the A(1905) resonance, the large 
difference between jjiN p and /^jvp can be attributed to a further suppression of the 
amplitude A(M^) due to the quark counting rules which require a 1/M® behavior of 
the RN^* vertex for the | ^ ^ transition. So, we have formulated the problem 
and outlined its possible solution. 

Note that interference between the different p-meson states does not change 
essentialy the dilepton contribution from the pion annihilation channel in heavy- 
ion collisions, since the p-meson form factor involved into that process falls off 
asymptotically like 1/t due to the quark counting rules, and so there should be no 
destructive interference between members of the p-meson family. The VMD model 
by KroU, Lee and Zumino allows two independent couplings for the photons and 
vector mesons and can be used to resolve the discrepancy between the photon and 
p-meson branchings of the nucleon resonances |36| . It provides, however, the form 
factors asymptotics F{t) = 0(1) in disagreement with the quark counting rules. 

3.2. Relativistic treatment of the vector meson decays of nucleon resonances 

The relativistic treatment of the R NV decays is in details discussed in [37]. 
Here, we sketch out the basic concept. For nucleon resonances with spin J > 1/2 
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and arbitrary parity, there exist three independent transition form factors, while for 
spin-1/2 resonances, two independent form factors should be considered 021E31^3- 
In terms of the electric (E), magnetic (M), and Coulomb (C) form factors, the 
decay widths of nucleon resonances with spin J = / + 1/2 into a virtual photon with 
mass M has the form j37[ I38j : 



I 



9a {l\f TOi(m|-M2)'+i/2(TO^_M2)'-i/2 



16 2'(2/ 

2 



G 



(±) 

M/E 



+ {l + l){l + 2) 



^E/M 



TO? ^ 



where to* refers to the nucleon resonance mass, to is the nucleon mass, m± = to, ±to. 
The signs ± refer to the natural parity (l/2~, 3/2+, 5/2~, ...) and abnormal parity 
(1/2+, 3/2^, 5/2+, ...) resonances. G^j^^ means G\.j or G^. The above equation is 
valid for / > 0. For / = ( J = 1/2), one gets 



G 



(±) 

E/M 



a 

2 



(to| 



Af2)3/2(™| -M2)1/2 



g: 



(±) 



(25) 



If the width T{N* N'y*) is known, the factorization prescription can be used 
to find the dilepton decay rate: 



dT{N* Ne+e~) = T{N* ~* Nj*)MT{j* 



where 



Mr(7* 



4to? 



(26) 



(27) 



is the decay width of a virtual photon 7* into the dilepton pair with invariant mass 
M. 

The couplings of the p- and uj- mesons and of their radial exitations can be taken 
into account to describe radiative and electroproduction helicity amplitudes satisfying 
the quark counting rules. The available data on the partial-wave analysis of the 
multichannel TrA^-scattering can also be included into the framework of the extended 
VMD model. In ref.[3S], the parameters of the model are fixed by fitting those data 
and by taking into account quark model predictions when the experimental data are 
not available. In this work, we use the model to calculate the dilepton production 
from the nucleon resonance decays. 

In the relativistic case, we cannot write a compact expression for the suppression 
factor (23). However, it is clear that the effect does exist. The destructive interference 
is prescribed by the quark counting rules which are implemented into the relativistic 
model. It means that the radiative decay should be less probable as compared to the 
naive VMD estimate from the ground-state p-meson. 

The medium can destroy the destructive interference, in which case one can expect 
an enhancement of the dilepton production below the p-meson peak [35] . A detailed 
treatment of this effect will be given elsewhere j39| . 
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The A(1232) resonance is treated in the same way as the other resonances. We 
take into acocunt 10 resonances hsted in Table 4. The p- and cj-meson channels 
are treated on the same footing. The numerical results demonstrate that besides 
the A^(1520) and A(1232) resonances, the iV(1535) and A(1620) have considerable 
contributions. It can be seen from Fig. 3 where the resonance contributions are shown 
for proton kinetic energy T = 1.61 GeV. At moderate invariant masses M < 0.35 GeV 
of the dilepton pair, the resonance contributions are dominated by the A(1232). At 
larger masses M > 0.35 GeV, contributions from the heavier resonances become 
dominant. 

4. Numerical results 

The results for the dilepton spectra are shown in Fig. 0] We show also inclusive and 
subthreshold cross sections separately. To compare with the experimental data, the 
acceptance of the DLS detector with respect to the e~^e~ pairs that have invariant 
mass M , transverse momentum p^, and rapidity y is taken into account. For each 
process, the distribution over the pt and y is determined by the available phase space 
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of the process and then weighted with the filter function f{M,pT,y) provided by the 
DLS collaboration. The details of this procedure are described in Appendix B. Finally, 
the finite mass resolution of the detector, AJW^^p* = ± 25 MeV, is taken into account 
by smearing the spectra with a Gaussian distribution which corresponds to a standard 
error of ct = 25 MeV. 

At the lowest initial kinetic energy of the proton, i.e. T = 1.04 GeV, the 
cross section is dominated by the 7r°-Dalitz decay below M < 100 MeV and by 
the Dalitz decays of the nucleon resonances, mainly the A(1232)-resonance, at 
M ^ 200 500 MeV. Compared to our calculation there is an excess of detected 
e+e~ pairs at M > 300 MeV. Earlier calculations |2S] obtained higher cross sections 
in this mass range. This is due to the normalization to the R N p branching ratios 
within the framework of the naive VMD which overestimates the radiative decay rates 
R — » A''7, as discussed in Sect. 3. 

At higher initial proton energies, the agreement with the DLS data is generally 
very reasonable. The contribution from the 77-meson Dalitz decay is dominant at 
M « 0.2 — 0.4 GeV, while the Dalitz decays of the baryon resonances dominate at 
M > 400 MeV. 

For the proton kinetic energy T — 2.09 GeV, the inclusive and subthreshold 
productions of the p- and w-mesons become important at AI w 800 MeV. The w-meson 
peak in the inclusive cross section is rather pronounced, whereas the subthreshold uj- 
meson production cross section is large and smooth. The total cross section, as a 
result, exceeds significantly the experimental data at M « 0.7 — 0.8 GeV. 

R,efs. |23l I25j agree quite well with the experimental data at the w-meson region. 
In ref.|2ni) the inclusive w-peak is not reproduced, apparently, due to a stronger ad hoc 
smearing with a mass- dependent parameter a = O.IM. The smearing at the w-meson 
peak turns out to be factor of 3 greater than it should {a = 80 MeV instead of 25 
MeV). In refs. |23[l^ . the w-chanels from the nucleon resonance decays are neglected. 
This is the reason for the difference between our results and relults of refs. 23, 25 . In 
our calculations, the w-chanel contribution is large. It can only be lowered by price 
of an ad hoc reduction of the decay probabilities of the nucleon resonances into the 
vector mesons, which has been calculated using the quark models 02] • The quark 
models reproduce, however, the p-meson decays of the nucleon resonances quite well, 
i.e. they are in good agreement with the partial-wave analysis of the ttN inelastic 
scattering 0^1 ■ In this way the experimental underestimation of the dilepton yield at 
the uj-peak at T = 2.09 GeV seems to contradict the irN inelastic scattering data. 

Remarkably, the three highest experimental points at T = 2.09 GeV lie above the 
kinematical limit Mmax = \/s — 2mN ~ 850 MeV and the indicated experimental error 
/^Mexpt = ± 25 MeV is not sufficient to explain their occurrence. We suppose that 
the experimental resolution is not good enough to resolve the kinematic threshold. 

Finally, at T = 4.88 GeV, the contribution from the inclusive production of 
the 77-, jO-, and w-mesons becomes dominant at M « 300 800 MeV. There is an 
underestimation of the dilepton yield in the region M k, 400 -f- 700 MeV. A similar 
underestimation was found in both, at T 2.09 GeV and T 4.88 GeV. As 
proposed in ref.[5S], the existing gap might be filled by the subthreshold dilepton 
production via the baryon resonances. However, we were not able to match the 
data using a consistent description of the photoproduction data and the R Np 
meson decay branchings, with the proper application of the Breit-Wigner formula 
(see details in Appendix 1), and removing possible sources for the double counting. 
Each of this three aspects leads to a reduction of the dilepton yield. Therefore, one 
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Figure 4. The differential dilepton production cross sections as a function of 
tfie dilepton invariant mass, M, after applying the experimental filter and the 
smearing procedure (see text). The solid curves are the total cross sections, the 
dashed curves correspond to the inclusive production, and the dotted curves 
correspond to the subthreshold production. The experimental data are from 

ref.irn. 



cannot exclude that the origin of the so caUed "DLS puzzle" can be traced back to 
the elementary pp level and is not a specific feature of heavy-ion collisions. New 
experimental measurements of the dilepton cross section, especially at T = 1.04, 2.09, 
and 4.88 GeV, would certainly help to clarify this point. 

In this context one should be aware that the comparison to the experimental 
data is strongly influenced by the acceptance of the DLS detector. In Appendix 2, we 
discuss the application of the corresponding filter program |221 for the calculation of 
the experimentally measured cross sections. In Fig|Slthe effective detector efficiency, 
smeared by a Gaussian distribution with the standard deviation cr = 25 MeV, is 
shown as a function of the dilepton mass M for decays tt 7e+e^, -q —t 7e+e~, and 
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p°(a;) e+e~ at the two highest proton energies T = 2.09 GeV and T = 4.88 GeV, 
where the effects of the multiple pion production are most important. It can be seen 
that the effective acceptance decreases with increasing energy for a fixed number of 
pions in the final state. On the other hand, when the number of pions increases, the 
acceptance increases as well. This can be interpreted to mean that a larger number of 
the pions reduces the available phase space for mesons decaying to the dilepton pairs, 
and the decays of such mesons can be detected with better efficiency. 

The effect is particularly strong for the tt '~fe~^e~ decay at T = 4.88 GeV. While 
the average pion multiplicity is around 2, the effective acceptance is extremely small 
below n « 6. The acceptance is not reliable when it is much smaller than unity. In our 
case this happens at n < 6. While the statistical distribution gives here, as we expect, 
reasonable estimates, the calculation of the part of the cross section connected to the 
additional production of n < 6 pions turns out to be unreliable. From the other side, 
at n > 6 the filter is well defined, but the binomial distribution gives exponentially 
small probabilities. The highest part of the pion spectra corresponding to n > 6 also 
cannot be calculated accurately. So, we consider the difference between our results for 
the pion contributions and those of refs.[2SII^ by about a factor of 3 and those of 
ref.|2ni by about a factor of 6 as a conservative estimate for uncertainties inherent in 
the theoretical calculations for both, the distribution over the pion multiplicities and 
the experimental filter acceptance in the region of small invariant masses. 

At lower energies, these uncertainties practically disappear. In FigEJ we show 
a plot for the tt 7e+e~ decay at T = 2.09 GeV. It can be seen that now already 
for n = the effective acceptance is no more extremely small. For heavier mesons, 
the calculation of the acceptance is safer, which is again connected with less energy 
available for the produced mesons and, respectively, a better efficiency for the detection 
of the dileptons. 

For the application of the filter we assumed an isotropic distribution of the 
particles in the final states in the cm. frame. This is justified at small energies 
T. With increasing kinetic energy, the distribution acquires a bias towards the beam 
direction. This is an additional source of uncertainties in the calculations of the filter, 
which can be important at energy T — 4.88 GeV for the pion Dalitz decays. 

The many-body phase spaces entering into Eq.(4) are known to be very sharp 
functions of the arguments. The Breit-Wigner distribution over the dilepton mass, M, 
gets therefore an enhancement towards small values of M. The greater the number of 
pions in the final state, the more important is this effect. We found that rare processes 
with probabilities w„ < 0.03 corresponding to large numbers of pions produced in 
association with the vector mesons, that should in principle give small contributions 
to the total cross sections, become very important at masses M < 200 300 MeV. 
The spectral functions of the vector mesons are not known well far away from the 
vector meson peaks. This effect is thus beyond the scope of the present model. It 
should be analysed separately. In the present calculations, we apply a 3% criterion 
to the multiple pion processes: The values Wn are set equal to zero every time when 
Wn < 0.03. This works for the inclusive vector meson production at T = 4.88 GeV. At 
smaller energies pions are not produced in association with vector mesons. 

5. Conclusion 

We have considered the dilepton production in pp collisions at BEVALAC energies. 
The subthreshold production of vector mesons through the nucleon resonances is 
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Figure 5. The effective acceptance of the DLS detector versus the invariant 
dilepton mass, M, for different numbers of pions, n, produced in association with 
the pion Tr'^, ry-meson, and p'^(aj)-mesons at two highest energies T = 2.09 and 
4.88 GeV. The numbers over the curves show the numbers n of the pions. 
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described within the extended VMD model which aUows to bring the transition form 
factors in agreement with the quark counting rules and provides an unified description 
of the photo- and electroproduction data, j{'y*)N N* , the vector meson decays, 
A^* Np{oj), and the dilepton decays, N* N£'^i~'. The dilepton decay rates are 
described relativistically using kinematically complete phenomenological expressions 
and numerical results of ref. [HHI- In this context, we discussed also the problem of 
double counting and proposed its possible solution. 

The resulting dilepton spectra are reasonably well described at proton energies 
ranging from T — 1.27-^ 1.85GeV. At T = 1.04 GeV, there exists an overestimation 
of the dilepton yield, at T = 2.09 GeV, we see an underestimation in the vicinity of 
the w-meson peak. At T = 4.88 GeV we observe an underestimation in the region 
of dilepton masses below the p-peak {M « 400 700 MeV). We hope that future 
experimental investigations will clarify these problems. 
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Appendix A. Breit-Wigner description of resonances v^rith 
energy-dependent widths 

Let us consider a process with a resonance in the intermediate state. It can be 
produced either in a two-body collision or as a result of the decay of a particle or 
another resonance. Let the resonance further decay to some specific channel i. The 
amplitude for the total process, M^, i.e. the amplitude for resonance production, 
propagation, and subsequent decay to the channel i is a product of the amplitude 
of its production Mp, the resonance propagator, and the amplitude of the resonance 
decay 

1 

' p'^ — rn^ + ' 

where p is the momentum of the resonance, m is its pole mass, S(p^) is the resonance 
self energy. The pole mass m is defined such that i?eS(TO^) = 0. In general, ReYiijP') 
starts with terms of the order 0{{p^ — m^Y). These terms are further neglected. The 
imaginary part of is equal to 

/mE(p2)^i^|Alj^|2ci,j^(p2) ^ (A.2) 

i 

where is the phase space for the resonance decay into a channel i. We can 

therefore write either 

/ml](p2) = [1 ^ I j^|'<f ^,(/) ) ^ Vp^tUp") (A.3) 
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or 



/mS(p2) = (^_L ^ |A^j,p$j,(p2) j ^ rnff„,(p2). (A.4) 

Both definitions of the total width, r^((p^) and r^((p^), can be used in the relativistic 
Breit-Wigner formula, but the width should be multiplied by the proper resonance 
masses, V?? (running mass) and m (pole mass), respectively. || The square of the 
amplitude M*, integrated over the phase space of the final particles with momenta, 
Pp and p^, and normalized by the corresponding factors for the initial particles, gives 
either a cross section (two initial particles) or a width (one initial particle). 
For the scattering problem, the cross section has the form 

where j is the flux of the incoming particles. 

Further, we introduce two (5-functions corresponding to momentum conservation 
in the processes of production and decay of the resonance and the running mass M of 
the intermediate resonance 

= / (pi + p2 _ Sfep^- - p)^ (4) _ j:ipi^)d^p J 5{KP - p^)dM^ ^"^-"^ 



and obtain 



where 



^ V(M2-m2)2 + (Mr«t(M))2 ^ • ' 



Mr,(M) = mf,(M) = i|A^^|'$^d(M2). (A.8) 

Here, one should also use the partial widths Ti{M), TiiM) with the proper masses 
M, m. The similar arguments apply in case of the decay problem. 



Appendix B. Effective filter function 

A comparison to the DLS data requires to take the experimental detector efficiency 
into account. For this purpose a filter function is provided by the DLS collaboration. 
In particular at large T (e.g. T = 4.88 GeV) this filter function is not a small correction 
to the theoretical calculations but is crucial for the comparison to data. Thus, in this 
appendix we discuss the influence of the detector filter in our analysis. 

In terms of the cm. frame momentum variables, the filter function can be 
rewritten as 

fiPT,y,M) = fip*T,y*+yc,M) (B.l) 

II In ref . 1241 the combination mpF''^ - (p^) (physical p-mcson mass and F without " ") has been 
substituted into the Breit-Wigner formula. Such a combination leads to an additional factor mp/M 
in the dilepton production cross section and, consequently, to an overestimation of the dilepton yield 
below the p-meson peak. We have brought the attention of the authors of ref. 1241 to this circumstance. 
Recently, a new paper appeared, ref. 1251 . where that inconsistency has been removed. 
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where = pt is the transverse momentum of the dilepton pair, is the rapidity of 
the cm. frame L* with respect to the laboratory frame L of the colliding nucleons, 



v. = (B.2) 

T is the proton kinetic energy in the L frame. The distribution of dilcptons in the 
cm. frame L* is isotropic. This is a universal feature which does not depend on the 
specific type of the reactions and is connected to the form of the cross section (4) only. 
So, we work with the filter function averaged over the angles in the L* frame: 

f{p*,y,,M)= f{p*^,y*+y^,M) (B.3) 



1 



where 



Pt — P stnt), 

y* = Lln( ^), 

" 2 e*-»r, 



and 



— p*cos'd, 
e* = ^jAP+p*"^. 

The problem reduces to finding the dilepton distribution over the dilepton momentum 
p* in the cm. frame L* . 

The probability distribution of the dilepton momentum in the L* frame for the 
direct decays V — > e+e~ is given by 

dW{p*) = ^«;„7^„$2(y^,M,Mjf)dMi$2+„(Mjf...) (B.4) 

n=0 

where 

The effective filter function can be calculated as follows 

rff{T,M) = Y,w^rn^f{T,M) = / dW{p*)f{p\y,,M). (B.5) 

The value of Mx is integrated out within the limits 2mM + n^T^ < Mx < y/s — M. 
The momentum of the dilepton pair p* = p*{yG, M, Mx) is given by 

P (Vs,A/,Mx) = T-p^ . 

2y/S 

A 100% detector efficiency would yield /(pr, y,M) = 1 and /"^^{T, M) = 1 in virtue 
of Eq.@. The values f^^^{T,M) are plotted in FiglHlfor the p°(w) ^ e+e~ decays 
for different values of n at energies T — 2.09 and 4.88 GeV. 

For the Dalitz decays A4 — *■ A4'e^e~ , the probability distribution for the dilepton 
energy e* in the cm. frame L* can be written as follows 

dW{e*) = / ^^^i / ^^^(e* - Vk)de* 

n=0 J i2mM+nfj.^)^ J -1 ^ 

X $2(V^,M>l,Mx)$2+„(Mx...). (B.6) 
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Here, k is the four momentum of the dilepton pair, /i^ is the mass of the meson M. in 
the reaction pp — > MX with the subsequent decay M — > A^'e+e". In the L* frame, 
the vector is defined as 77 = (1, 0). The coefficients _D„ are defined as before. 

Now we should pass to the rest frame L** of the meson M , where 77 = 7(1, — nu). 
The 7-factor and the velocity v of the meson Ai in the L* frame are determined by 
equations 

7M^ ={s + Ai- M|-)/(2^/i), 

The unit vector n shows in the direction of the meson velocity in the L* frame. In 
the meson rest frame, L**, the dilepton pair has momentum k = (e**,n'p**) where 

p** =p*{fij^,^j^,,M). 

The function p* (...) is defined earlier. The unit vector n' shows in the direction of the 
dilepton pair momentum in the L** frame. In Ea. ljB.6| the value x is the angle between 
the directions of the meson velocity in the L* frame and velocity of the dilepton pair 
in the L** frame, so that cosx = nn' and therefore rjk = j{e** ,vp**cosx)- 

In Eq.(40) the integral over the angle x is evaluated explicitly, and we obtain 

= J2n^„D^-P^ / dMj,9ie\Mx)<P2+n{Mx-).iB.7) 

where 

6(t* Mx) I ^' ^^^** ~ - ^^'^** + ^P**^' (B 8) 

1 0; otherwise. 



The effective filter function can now be calculated to be 

r"{T,M) ^yw„f^"{T,M) = / dW{e*)fip*,y,,M). (B.9) 



n=0 



The values f^^^ {T, M) are plotted in Fig [51 for the 7r'^(ry) ^ e'^e~7 decays for 
different values of n at energies T = 2.09 and 4.88 GeV. 

It is now sufficient to multiply the differential cross section Q with the 
corresponding effective filter function f^-^f{T,M) < 1 in order to compare the 
calculations with the experiment. For the evaluation of the direct contributions one 
should use expression (jB.5p . while for the Dalitz decays one should use expression 
l)B.9|l . The function f^-^^ {T, M) is given by a two-dimensional integral in Ea. (|B7 
and by a three-dimensional integral in Ea. (|B.9|l . 
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